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Abstract 

Using a local moment approach of Logan et al. we developed a solver for a multi-orbital single impurity Ander- 
son model. The existence of the local moments is taken from the outset and their values are determined through 
variational principle by minimizing the corresponding ground state energy. The method is used to solve the dynam- 
ical mean-field equations for the multi-orbital Hubbard model. In particular, the Mott-Hubbard metal-insulator 
transition is addressed within this approach. 
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The problem of a magnetic impurity in a metallic 
host, i.e. the Kondo problem, plays a major role in 
many branches of modern condensed matter physics. 
This problem of coupling degenerate spin degrees of 
freedom to gapless fermionic excitations is essential in 
the dynamical mean-field theory (DMFT) of strongly 
correlated electron systems [1]. The DMFT is a self- 
consistent method, where in each iteration loop the 
Kondo problem is solved numerically. The application 
of the DMFT to investigate real systems, in particular 
with many orbitals, is difficult because of the neces- 
sity to solve the Kondo problem for arbitrary parame- 
ters. Numerically exact methods, like numerical renor- 
malization group [2] or determinant quantum Monte 
Carlo [3] are very time (CPU) consuming, in particular 
when the number of orbitals is large. Reliable analyt- 
ical methods are therefore needed. One of such meth- 
ods, which recovers Luttinger property and Kondo ex- 
ponential scale as well as Hubbard satellites has been 
recently invented by Logan et al. [4] and named a local 
moment approach (LMA). In this contribution we gen- 



Corrcsponding author. 
Tcl:+48 (22) 5532315 fax: +48 (22) 6219 475 
Email address: anna.kauch@fuw.edu.pl (A. Kauch) 



eralize this method to investigate multi-orbital systems 
at zero temperature. Within the DMFT framework 
supplemented with the LMA we analyze the problem 
of Mott-Hubbard metal-insulator transition (MIT) in 
orbitally degenerate and non-degenerate systems. 

We start with the generalization of the LMA to 
solve the single impurity Anderson model (SI AM) 
with many orbital levels: 
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where the direct U and U' as well as exchange J inter- 
actions between electrons of spin a and on orbitals a 
or /3 are taken into account. 

In the unrestricted Hartree-Fock approximation to 
the SIAM the local Green functions GZ 13 {^ a , ^p;^) HF 
depend explicitly on the local magnetic moments \x a 
on each orbital. These local Green functions are used 
in the random phase approximation (RPA) to deter- 
mine the polarization diagrams n°f (fi a , np\ ui), cor- 
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responding to spin flip processes, and hence the self- 
energies S"' 3 (/i cl , fj,p; to). The Green functions calcu- 
lated in this way depend now on the local moments 
fi a . However, in a single impurity problem the spin- 
rotational symmetry cannot be spontaneously broken. 
To restore this symmetry within the approximate so- 
lution we use symmetrized Green functions GZ^ (u) = 
2l^E±| Ma |,±| M(a | G^ific, fin; uj), where ±\fi a \ in the 
sum means summation over different directions of lo- 
cal orbital moments and N or t, is the number of orbitals. 
In this method the existence of the local moments is 
assumed from the beginning. 

The lengths of the local moments \/j, a | and the num- 
ber of particles on each orbital n a are determined 
by minimizing the ground state energy function, i.e. 
-^physical = mm {Mc«,«c,} E(fj, a ,n a ). In this work we also 
adjust the Fermi energy of the electrons in the impu- 
rity atom to have the proper total number of electrons. 

This generalized LMA is used to solve multi-orbital 
Hubbard model if H ubb = Y, a ,v ^tjAa^i^ + 
J^iocai, where the local part is a sum over all lattice site 
terms which are of the same form as the atomic part 
in the SIAM. This model is solved within the DMFT 
where the self-consistency condition relates the local 
matrix Green functions with the matrix of the self- 
energies [3]. This condition simplifies for the Bethe 
lattice that is used in this contribution. The band 
width W a = 2 for each orbital a sets the energy units. 

In Fig. 1 the spectral functions for two-orbital Hub- 
bard model are shown for degenerate case 6a = es (up- 
per panel) and for nondegenerate (asymmetric) case 
\ca — es| = 0.3 (lower panel). Also the cases with (J = 
J7/4) and without (J = 0) Hund exchange coupling are 
compared. The inter-band interaction U' = U — 2J is 
fixed preserving £?f/(4) symmetry. In all cases the av- 
erage number of electrons is fixed to n e = 2 per lattice 
site. 

In the first case = es when U is sufficiently large 
the correlation gap is opened and the system is a Mott 
insulator. At smaller U, when the system is on the 
metallic side, the Kondo-like peak is seen at J = and 
is destroyed at J > 0. When the orbital degeneracy 
is removed, as in \sa ~ es| = 0.3 case, the system for 
a given U is metallic on one of the orbitals, whereas 
on the other it is insulating. The number of electrons 
per site on each orbital level is no longer equal and 
the spectral functions on different orbitals begin to dif- 
fer qualitatively. At larger U the system is insulating 
in both orbitals. So we find an orbital selective Mott- 
Hubbard MIT in asymmetric case. Increasing further 
the difference \ea — es| would lead to the appearance 
of the band insulator possessing a gap between bands. 
Initially the orbital selective Mott-Hubbard MIT was 
found in a two-orbital symmetric case, where orbitals 
had the same energies but different bandwidths [5,6,7]. 
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Fig. 1. Spectral functions for two-orbital Hubbard model on 
Bcthc lattice with infinite coordination number. Upper panel: 
degenerate case with U — 2.0; solid line is for J — 0, dashed 
line is for J — t//4. Lower panel: nondegenerate case with 
\eA — cb — 0.3, U — 1.8 and J — U / '4; solid and dashed lines 
correspond to different orbitals. In both cases the Fermi level 
is at zero energy 

The present contribution extends this investigation on 
a case where this orbital degeneracy is removed. The 
crystal level splitting can be tuned in some extend by 
the axial pressure and thereby there is a practical op- 
portunity to search experimentally for the orbital se- 
lective MIT. 

The generalized LMA to multi-orbital SIAM allows 
us to efficiently study the correlated electron systems 
within DMFT. In particular it is relatively easy to ad- 
dress the problems where the degeneracy of orbital 
levels is removed, and hence the Fermi energy has to 
be determined self-consistently in each iteration loop. 
Also the local moment approach is much more efficient 
in studying problems with more then two orbitals [8] . 
This regime is hardly accessible in numerically exact 
approaches. 
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